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Dipole -^change  Soin  Waves  in  a Cylindrical 
Ferromagnet 

Lai  Wu-yan  Wang  Ding-sheng  Pn  Pn-cho 

(institute  of  Physics,  Academia  Sinica) 

Abstract 

Taking  both  dipole  and  exchange  interaction  into  account, we  have 
deduced  the  secular  equation  for  the  dipole-exchange  spin  waves  in  an 
axiall,r  magnetized  cylindrical  ferroma'met.  Numerical  results  o^  frequencv 
spectra  are  also  given.  When  the  wave  vector  the  spin  wave  is  small, dipole 
energy  is  the  most  important  term  and  our  results  coincide  with  that  of 
magnetostatic  modes.  When  the  wave  vector  becomes  lar^e,  exchange  energy 
plays  a great  role,  and  our  results  approach  those  of  the  theory  of 
marcroscopic  exchange  waves. 


1.  Introduction 

This  article  is  intended  to  discuss  the  behavior  of  spin  waves  in 
an  axially  magnetized  cylindrical  ferromagnet  when  the  dipole  and  exchange 
energy  is  of  approximate  magnitude.  We  call  this  kind  of  waves  as 
dipole-exchange  spin  waves.  For  instance,  in  yttrium  iron  garbet  (YIG) , 

n *i 

when  wave  vector  Bja£- 10  cm  , these  two  sets  of  energy  are  almost  the  same. 

In  a region  where  lO^crrf 1 , diploe  energy  is  the  major  one;  it 

is  the  so-called  magnetostatic  mode  rerion.  Joseph  and  Schlomann  have 
defined  the  magnetostatic  mode  in  cylimtrical  ferroma'met.  In  a region 
where  IO^cit T-,  exchange  energy  is  the  major  one*  this  is  the  so-called 

exchange  mode  region.  Herring  and  Kittel  have  found  the  spin  wave 
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frequency  spectra  by  using  semi-classical  met’uod  in  their  study  of  the 
theory  of  macroscopic  spin  waves.  It  can  be  possible  to  discuss  frequency 
sp  'Ctra  within  the  range  of  a whole  wave  length  only  if  both  dipole  and 
exchange  interaction  is  taken  into  account  at  the  same  time.  The  solution 
o'*  energy  spectrum  of  dipole-exchange  spin  waves  of  plane  r‘iln  fcrroma~nct 
can  be  found  only  in  de  Wames  and  Wolfram's  works. 

Through  discussion  of  dipole-exchange  spin  waves  in  a cyclindrical 
f erroma'Tict,  in  this  article,  we  can  find  frequency  snectra  of  whole  spin 
wave,  '//hen  spin  wave  vector  Y is  smaller,  it  is  corres  ondent  to 
magnetostatic  node.  When  is  larger,  it  approaches  the  result  of  exchange 
mode  theory  of  macroscopic  spin  waves#  In  the  middle  region,  spin  wave- 
is  a mixture  of  magnetostatic  mode  and  exchange  mode,  and  at  the  same  time 
it  has  body  wave  and  surface  wave.  State  density  has  noticeable  changes, 
it  cannot  simplv  add  magnetostatic  ener-y  and  exchange  energy  up  to 
substitute  for  whole  freouency  spectra. 


2.  Basic  Methods 

Here  is  a dipole-exchange  spin  wave  model  in  an  ax'allv  magnetized 
cylinder.  When  the  wave  length  of  spin  waves  is  much  louver  than  atomic 
separation,  the  sp'n  waves  can  be  described  by  motion  equation  of  macroscopic 
magnetic  moment.  And  the  motion  of  macroscopic  magnetic  moment  can  satisfy 
magnetostatic  equation  (omit  dissemination  effect) 


>-o, 

+ divm  — 0, 

and  motion  aquation  of  magnetic  moment 

±M--rMxH . 


(2.1) 

(2.2) 


In  it,  effective  field  H~Ht  + h + — DViM  is  the  sum  of  constant 

A*» 

megnetic  field  FQ  (along  cylindrical  axis  z),  microscopic  magnetic  filed  h 

equavalont  exchange  field  “ DV’M  (d  is  the  parameter  of  exchange 

h* 

action,  for  example,  in  YIG,  D = 2.6  x 1 0~^cm'’). 

The  relation  between  magnetic  moment  and  time  is  e’m  , taking  t ' me 
factor  off,  equation  (2.3)  is  then  written  in  component  form  as 

w,0  ~ <*A)mx  — jtom , “ fhum^>r  (2.4) 

jwm,  + c0o(l  — aA)m,  “ (2.5) 

in  it  "0  “ rW,  a»«  — rM,/fh,  * — Dum/u» . 

Adding  scalar  quantity  magnetic  potential  <A,  and  male  in" 
h — — V</»,  then  from  equations  (2.l)-(2.3),  under  linear  approximation,  a 
sextic  differential  equation  of  the  internal  magnetic  potential'/*1’1  of  the 
sample  can  be  produced 

[(£>  - & + e)v>  - e - o,  (2.6) 

in  it©  is  a differential  operator,  0 — — aV’),  and  con8'tant 

Q — u/um,  Qm  » 

On  the  outside  of  the  sample,  magnetic  potential  <pw  satisfies 

eluation  - o.  (2  7) 

For  a axially  magnetized  cylinder,  the  solution  of  equation  (2.6) 
should  be  of  the  following  form: 

<*"  - «p). 


(2.8 


Substituting  equation  (2.2)  for  equation  (2.6),  we  have  <Pt(p,  <p),  by 


which  the  equation  can  be  satisfied  is 

(A,  + )(A,  + 4i)(A,  + $ )0f(p,  V ) “ 0,  (2.9) 

in  it  A<  ” T 0*  -J-)  + “3  & & -tf  are  three  roots  of  a cubic 

P Op  \ Op/  p1  Op* 

algebra  equation  of  k 

(*  + py  + ±+Q*  (*  + py  + LzJPOL  + Om-  *FQm  ^ + ^ _ ££*  . 0 

a tP  B* 

(2.10) 


It  is  easy  to  see  that  the  solution  of  equation  (2.9)  is  the  sum 
of  solutions  of  three  'Telmholtz  equations  (A,  + $)#  «■  0,  (»  “ 1,2,3)  • 
thereupon  the  "eneral  solution  of  equation  (2.6)  is 

* • - S ("  ^ S p<*>  cm) 

• 4-1 

in  it  a is  radius  of  the  cylinder  and  J,(.kiP')  is " Wt  (»  ""  0,  ±1,  ±2-*»)  Bessel 


function. 

solution, 


The  items  indicated  by  Neumann  function  Nn(klp)  ^Ln  general 
because  of  the  limited  boundary  condition  <pl*(.p  “■  0^  are  cxluded. 


Fbr  the  outside  of  the  sample,  the  general  solution  of  equation 


(2.2)  is 


**"  - 2 ("  W-KMW—"*  (2.12) 

in  it  kXPp)  is  transformed  Bessel  function  of  the  second  kind.  /.(fr>) 
indicated  by  transformed  Bessel  function  of  first  kind  in  general  solution, 
because  boundary  condition  </’"‘,(p_*  00 ) "“  0 , are  excluded. 


3.  Boundary  Condition 

The  corresponding  boundary  conditions  a “ * of  magnetostatic 
equations  (2.1)  and  (2.2)  are  tangential  component  succession  of  h and 
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normal  component  succession  of  <“o h + nt,  namely 


*toU-*“U,  (3.!) 

(Mi"  + ^oAr | (3  2) 

For  nameto  static  mode,  <*  ■=  o,  equations  (2.4)  and  (2.5)  are  transformed 
into  alaebra  equation?,  and  only  equations  (2.1)  and  (2.2)  remain  to  be 
differential  equations,  and  there  are  only  two  coefficients,  in  their 
corresponding  solutions.  The  mametostatic  secular  model  can  be  determined 
by  these  two  boundary  conditions. 

For  dipole- exchange  mode,  there  must  be  boundary  conditions 
correspond ' ng  to  differential  equations  (2.4)  and  (2.5).  There  is  a frequently 
used  boundary  condition  which  assumes  that  the  boundary  magnetic  moment  is 
same  as  the  internal  magnetic  moment,  and  can  satisfy  same  motion  equation. 

This  raea:'S  that  the  surface  magnetic  moment  is  completely  ".free”.  Thus 
from  equations  (2.4)  and  (2.5)  and  by  using  generalized  function  method, 
if  boundary  plane  is  proved  to  be  a cylindrical  plane,  the  "free1*  boundary 
condition  is 

Another  frequently  used  boundary  condition  makes  , ***(p  This  means 

that, because  of  some  reason,  it  makes  magnetic  moment  of  the  boundary 

plane  completely  "nailiiig  upon"z  direction.  It  may  use  a certain  combination 
(partial 

A nailing  up)  of  the  two  conditions.  In  Document  (3),  in  the  discussion  of 
dipole-exchange  spin  waves  of  plane  film,  the  surface  is  considered  as  a 
crystal  plane.  Then  through  the  process  of  microscopic  motion  equation 
approaching  to  succession,  a corresponding  boundary  condition  is  obtained. 
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The  merit  of  such  a solution  is  that  when  it  is  short  wave,  it  is  possible 
to  approach  an  exchange  mode  that  is  acquired  through  the  solution  of  a 
microscopic  motion  equation.  But  it  is  difficult  to  extend  this  method 
to  cylindrical  surface.  The  operation  that  follows  is  to  use  "free" 
boundary  condition,  and  +he  impact  of  boundary  conditions  will  be  discussed 
in  Section  5. 


4.  Secular  Equation 

The  solution  of  equations  (2.11)  and  (2.12)  can  satisfy  equation 
(2.6).  From  A — — V</»,  it  s possible  to  secure  a corresponding  h.  To 
seek  m from  h,  it  should  be  under  boundary  equation  (3.3)  to  try  to  solve 
differential  equations  (2.4)  and  (2.5).  Using  Green  function  method,  it 
can  have 


(I)) 2 

/ 1 -i\ 

x \-i  -t)  u-.ap)  - 

in  it 

i)-^,  n1-  <.QQm+ 

a * 

b*  _ *.«/.+ .g-  _ ki*f~-\(.kia')  + Jn-^ki*). 

+ ■/.+,(£«)  * l«)  + Jn-xiv'i 


In  order  to  make  solution  of  equation  (4.1)  satisfy  equation  (2.2) 
can  prove  hat  coefficient,  should  be  made  satisfy 


» 


it 


(4.2) 

(4J) 
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f 

X. 


(4.4) 


By  substituting  equations'  (2, ll) , (2.12)  and  (4.4)  for  boundary  condition 
equation  (3.1)  2nd  (3.2),  the  coefficient  of  & acquired  should  satisfy 
the  other  two  conditions  which  are 


2 - Dm,Km(Ja\  (4-5) 

i m\  . * . 


2 cM*Uki*) 
» » 1 


fiakMXA.S') 


The  nonrzero  solution  condition  of  the  linear  equations  (4.?),  (4.3), 

(4.5)  and  (4.6)  of  coefficient  of  ^ is  the  secular  equation  that  determines 

dipole-exchange  node . 


Prora  equations  (4»2)  and  (4.3)  it  can  have 


c\,  _ I(L_L2 . Mii£l . ci,, 

E(2,  3) 

d «.  . cl, 

E(  2,3)  ;.<*,«)  •** 


r 

Substituting  the  above  equation  for  (4.5)  and  (4.6),  it  can  have  linear 

equation  of  C'.,  HD.,,  The  non-zero  solution  condition  is  y 


in  it  _ v , ' 

EC../)  - | % g‘\’  . : 

f,  - »(*  - 0^4^  + ia-y  - ^ 

c,  - nXfidUkiLk-  - 1),  ■ 

Mi  “ —p1/  » 

X,  - DQm 0 - /*,)/[  1 + a(tf  + /?’)]  W-  1 »2»3-  . 

a*vL 

Inder  the  given  Qm,  a,  pJtfin  , Q which  can  satisfy  equation  (4.8)  is 
he  secular  frequency  of  dipole— exchange  mode. 


5.  Results  and  Discussions. ' 

The  three  roo!s  of  equation  (2.10)  corresponding  to  potential  equation 
(2.11)  and  mametic  moment  distribution  equation  (4.4)  have  three  oscillations, 
and  the  secular  equation  of  dipole-exchange  mode  is  the  accumulation  of 
these  three  oscillations.  When  Q<Q~'(.  1 + a/?2)  V ft  & ' kl  are  all 
negative  . But  within  tins  range,  to  secular  equation  (4.8)  there  is  no 

' , A 

solution.  Whan  O > 1 + ap1')  < 0 , and.  this  indicates  that 

of  the  three  oscillations  which  are  included  in  the  secular  equation,  one 
is  body  mode  and  the  other  two  are  surface  mode. 

Under  the  limit  of  long  wave,  when  p l/V  <*  > and  to  the  range  o' 

Q~'(l  +ap')<Q<  £K  (1  + Omyn  $ *•  fto*  - + 0?  - . This  moans 

v . > * 

that  k]  and  the  magnetostatic  body  mode  are  the  same.  To  the  range  of 
a >it'(i  + Qmy*t  ft  PK&  - XEO/CJE?  + 02  - G1)  . This  indicates  that  k2 
approaches  magnetostatic  surface  mode. 
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k-j  represents  one  surface  mode,  and  it  will  quickly  decay  following 
the  increase  of  its  d? stance  from  the  surface.  Because  $ < $(£ 

.-(^  + '1/- +■«'/-).  , it  concentrates  only  in  the  surface  layer,  of  which  the 
thickness  in  about  if ^Q;  level  of  quantity.  Among  the  magnetostatic  modes, 
there  is  no  one  which  can  be  same  as  this  one.  \J~&QC  has  a length  not,  which 
marks  the  characteristic  length  of  exchange  action.  To  YIG,  when  Qm—> 2t 
- jUj)  - 2.3  X 10-‘cm. 


Figure  1 is  the  relationship  of  which  are  secured  from 

the  solution  of  equation  (2.1 0),  when  £„~2,  fi  - 0.1(1  a/T)  . Figure  2 in  a 
contrast  between  i(  and  — + — Q 0 of 

magietostatic  mode  (the  symbols  used  in  Figure  1 and  Figure  2 are? 

They  clearly  represent  the  results  o^  the  kjre^oing  discussions. 


Figure  1 The  three  oscillations 
of  dipole-exchange 
mode 

0 — O.I(l/v'«’) 


Fi.gure  2 Comparison  of 

dipole-exhance  mode 
and  magnetostatic 
mode 

a. -2,  0-0.  !(!/✓«) 
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Under  riven  and  n,  the  solution  of  secular  equation  (4.0)  rives  a series 
of  eigenvalue . 0,„t  r-  1,  2,  3-..  is  the  order  number  o'*  eigenvalue,  and 
different  r rep?  eser.ts  different  zero  number  which  occurs  when  potential 
function  1/?  at  the  region  where  p •*  [0,«]. 


Firure  3 Dipole-exchange  mode 
frequency  spectra 
(solid  line) 

O.  - 2,  t - 50/a,  ■ « 1 


Figure  3 is  the  eigenvalue 

Ctm£r  ■»  1,2,*  • -7),  of  dipole-exchange 
mode  obtained  from  the  solution  of 


equation  (4.0)  when  Qm  — 2, 
a — 50\/ a , n ■=  1 . At  the  same  time, 
it  fives  spin  wave  (exchange  mode) 
frequency  spectra  of  macroscopic 
spin  theory.  For  the  situation  of 


wave  vector  paralleling  w!th  magnetic 


direction 


Q,  + of). 


for  comparison,  it  "ives  three 


curves  (r  — 1,2,3)  in  magnetostatic 


modes  of  cylindrical  f erromarmet . 


At  short  wave,  when  /?  n(l/\/  a) 

and  f(l/V  a),  the  theoretical 
results  of  dipole-exchange  mode  and 


' 

' 
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energy  spectra  and  magnetostatic  nodes  are  the  sane.  When  the  cylindrical 
radius  is  not  ver y large  (as  « — 50\/ o in  Fimire  3)>  tv,e  impact  of  exchange 
energy  can  already  be  seen.  And  t e ’ipole-exeMn *e  modes  are  probably  alread; 


in  shape  of 

— Q,„  (nn'metost-g  !c)  + exchange  correction  (5.1) 

p 

Ir.  a re  ion  where  l/a  « s « 1/v  o , k^  acquired  from  the  solution  of 
equation  (2.10)  includes  li?l>p*»i/a>  So  oscillation  kp  can  only  have 
some  action  on  the  surface  of  the  sample,  and  its  contribution  to  exchange 
energy  is  not  great.  The  exchange  correction  in  equation  (5.1)  comes 
mainly  from  k]_,  so  exchange  correction 


Roughly  sneakin',  t e distribution  of  eigenvalue  probably  is  + — 

so, only  if  * » a and  under  the  situation  that  n and  r are  not  vrey  large, 
the  correction  is  but  a small  one.  This  is  just  the  situation  indicated  in 
Figure  3. 


The  conventional  methods  are  simply  to  take  spin  wave  spectra  as  an 
accumulation  of  magnetostatic  energy  spectra  and  exchange  energy  spectra, 
and  rive  an  analome  result  to  equation  (5.1).  tv0m  strict  computation, 
it  is  known  that  it  can  be  established  only  when  is  larger  (short  wave) 

/t  / 

or  when  W is  smaller  (ion"  wave).  Then  t^ere  is  no  meat  difference  from 
the  rigid  theories.  But  this  cannot  be  applied  to  the  medium  range, 

P ~ 104 — 10 '/cm,  namely  dipole-exchange  mode  range. 


When  P < l/«,  the  difference  between  dipole-exchaneg  mode  and 
magnetostatic  mode  is  very  clear.  In  magnetostatic  mode,  when  p-+0t  all 


frequency  spectra  of  n and  r will  be 


o,„  (tem.)  u - & (i  + amy*. 

After  considering  exchange  action,  there  will  be  no  more  combination  of 

different  n and  r.  If  the  cylindrical  radius  is  very  largo,  to  the  3mall 

n and  r,  bV.s  split  is  not  noticeable.  3ut  to  radius  a, which  is  not  large, 

as  a -»  50%/ « in  Firure  3,  this  phenomenon  cannot  be  overlooked.  Besides, 

in  t e range  of  P 1/a  , w1  en  Q . is  large  and  U!l  has  been  reduced  to 

below  l/a',  so  the  exchange  energy  of  surface  mode  kp  becomes  important. 

After  considering  exchange  action,  secular  function  is  always  a mixture  of 

body  mode  and  surface  mode  and  this  kind  of  mixture  will  become  especially 

of 

remarkable  when  intersecting  /-  body  mode  and  surface  mode  takes  place.  In 
, surface  mode  is  no  longer  a single  curve,  and  it  has  been  cut  into 
several  sections  which  connect  the  body  mode  (Figure  3).  Talcing  dipole- 
exchange  mode  r = 5 as  example,  in  section  A,  body  mode  (kj)  is  the  major 
part  nixed  with  surf-ice  mode;  in  section  B,  surface  mode  (kp)  is  the  major 
part  mixed  with  body  mode;  in  section  C,  body  mode  is  the  major  part,  and 

Jt 

following  the  increase  otjtf  the  surface  mode  is  continuously  reduced;  in 
section  D,  when  it  approached  short  wave,  exchange  energy,  after  surpassing 
magnetostatic  energy,  gradually  becomes  into  exchange  body  mode. 
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